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Abstract: A clinical trial on multi-spectral imaging of malignant and non- 
malignant skin pathologies comprising 17 melanomas and 65 pigmented 
common nevi was performed. Optical density data of skin pathologies were 
obtained in the spectral range 450-950 nm using the multispectral camera 
Nuance EX. An image parameter and maps capable of distinguishing 
melanoma from pigmented nevi were proposed. The diagnostic criterion is 
based on skin optical density differences at three fixed wavelengths: 540nm, 
650nm and 950nm. The sensitivity and specificity of this method were 
estimated to be 94% and 89%, respectively. The proposed methodology and 
potential clinical applications are discussed. 

© 2012 Optical Society of America 

OCIS codes: (110.4234) Multispectral and hyperspectral imaging; (170.4580) Optical 
diagnostics for medicine. 
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1. Introduction 

Cutaneous melanoma currently represents 5% of newly diagnosed cancers in men and 6% in 
women; it is the leading fatal illness affecting the skin and is responsible for 80% of deaths 
from skin cancer [1]. It is hard to distinguish visually melanomas from non-malignant nevi. 
By now biopsy sampling from melanoma (which may increase the risk of spreading 
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metastases [2]) is the main method for making a precise diagnosis of the malformation. More 
patient-friendly non-invasive primary diagnostics would be preferable. 

Multispectral imaging (MSI) is a non-contact, optical method that can be used in vivo with 
potentially promising applications to skin studies and pathology diagnosis [3]. Using narrow- 
band filtering at different wavelengths, a set of spectral images of the skin/pathology is taken; 
spectral filtering is usually performed by Liquid Crystal Tunable Filters or a set of discrete 
band pass filters. Each pixel of the image contains information about the corresponding skin 
spot at a specific wavelength, and a diffuse reflection spectrum for a single pixel can be 
obtained by processing the set of multi-spectral images. The penetration depth of optical 
radiation in the skin tissues depends on wavelength [4], so diffuse reflectance from skin at 
broadband illumination may provide morphological information from different depths, and 
multispectral imaging analysis makes it possible to extract this information [3,4]. 

The optical properties of melanin and of the different stages in a potentially malignant 
histological transition (dysplastic nevi — » melanoma in situ — » melanoma) have been 
investigated by contact diffuse reflectance [5]. Optical properties of benign and malignant 
pigmented lesions have also been measured by oblique incidence diffuse reflectance 
spectrometry [6]. Multispectral imaging provides spectral information over the entire area of 
the lesion and has potential for in vivo skin diagnostics, including melanoma diagnostics [7]. 
Previous investigations indicate that reflectance of the benign nevus at 940 nm is close to that 
of the surrounding skin, whereas pigmentation of the melanoma is detectable [3]. 
Consequently, combined visible and near infrared image analysis might help to distinguish 
melanomas from other skin malformations. In this study multispectral images were taken in 
the spectral range 450-950 nm in order to find a suitable algorithm that could be used for 
melanoma-nevus differentiation. 

2. Methods and equipment 

2. 1 Experimental setup and processing algorithm 

The measurement setup is shown in Fig. 1. The skin area to be studied is illuminated by a ring 
of halogen lamps. The light from the halogen lamps passes through a diffuser and a polarizing 
film. Reflected and back- scattered light from the skin passes through an objective lens into a 
multispectral imaging camera (Nuance EX). The Nuance Imaging Module contained the 
principal imaging components in a single compact enclosure: a high-resolution, scientific- 
grade CCD imaging sensor, a solid-state liquid crystal filter with a polarizer, a wavelength 
tuning element, a spectrally optimized lens and internal optics [8]. To avoid artifacts caused 
by light reflection from the skin surface, the polarization of the illuminating light was 
orthogonal to the orientation of polarizer located in the Nuance EX camera. Besides, the 




- Polarizer 



Fig. 1. The measurement setup: multispectral camera Nuance EX, halogen light source, 
polarizer, and computer. 



#158653 - $15.00 USD Received 23 Nov 201 1; revised 18 Jan 2012; accepted 18 Jan 2012; published 9 Feb 2012 
(C) 2012 OSA 1 March 2012 / Vol 3, No. 3 / BIOMEDICAL OPTICS EXPRESS 468 



halogen lamps were powered from a stabilized direct current source to avoid artifacts due to 
temporal fluctuations of the illumination. The Nuance program CRi was used to acquire 
image cubes and to average the spectra obtained from pigmented and normal skin areas. 

During a clinical trial 82 pigmented lesions were studied by multispectral imaging in the 
wavelength range 450-950 nm. Spectra were obtained at different wavelengths over this range 
in steps of 10 nm. The half-b and width of the color filter that regulated each step was 15 nm. 
Overall, 17 histologically confirmed melanomas and 65 nevi identified by dermatologists 
were investigated in this study. The melanoma had different thickness: 8 melanomas fell in 
the thickness range from 1 to 2 mm and 9 within the thickness range of 3-5 mm (Breslow's 
depth). The nevi had diameters of at least 5 mm. 51 images at different wavelength bands 
were taken from each pigmented malformation. Before the measurement, the lesion was 
covered by white paper and a multispectral image set was taken, in order to record the 
illumination spectrum. Since white paper is not an ideal diffuse reflector, the spectral 
reflectance value was corrected by comparing it with an Avantes WS-2 reference tile (99%).. 
The respective corrections were further added to the captured image array using the same 
reference units (r.u.) that were also used for processing the pathology images. 

All the measurements were taken from patients in supine position, in order to minimize the 
motion artifacts; in addition, the image stabilization program Image J [9] with plugin 
TurboReg [10] was used. During the image acquisition process the patients were asked to hold 
their breath about 30 seconds while the image cube was taken. The duration of the whole 
procedure depended on the location of malformation. The following steps of the procedure 
were taken: 

• Patient relaxing in a comfortable recumbent position (about 1-3 minutes); 

• Acquisition of the reference and image cube (about 1-2 minutes); 

• Stabilization of the image cube, if necessary (about 5 minutes); 

• Image processing (about 15 seconds). 

The entire diagnostic process took about 15-20 minutes. 

The optical density was calculated by the CRi Nuance program as 

OD(A) = -log(I(A)/I 0 (A)\ (1) 

where l(X) represents the intensity of the skin-reflected light and lo(X), the intensity of the light 
reflected from the white etalon (at the same distance); X denotes the wavelength. 

In order to differentiate between melanoma and nevus, the following differentiation 
parameter p was proposed: 



p = OD 650 +OD 950 -OD 540 =\og 

where 



k-- 



7(540) 



7(650) -7(950) 



(2) 



fc = Z 0 (650)-Z 0 (950) 
/ 0 (540) 

is the coefficient that describes the reference, and OD 540 , OD 650 and OD 950 are the optical 
densities at 540nm, 650nm and 950 nm, respectively. The wavelength 540 nm corresponds to 
the wavelength of maximum absorption of the blood hemoglobin [11], while at 650 nm 
hemoglobin absorption is minimal and melanin absorption is more pronounced [12], and 950 
nm is the wavelength at which light penetrates the skin surface most deeply [3] in our 
measurement conditions. The parameter p was calculated for each pixel, and then parametric 
/7-images were constructed. 
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Figure 2 presents the patient-averaged spectra of the optical density for healthy skin (dark 
filled circles), common nevi (filled circles) and all melanomas (empty circles). Error bars 
representing one standard deviation are presented at few points of the curves to illustrate the 
statistical variation between patients. Each point of the mean OD spectra was calculated as the 
average value of all pixels in manually selected rectangular areas of regions of interest within 
and outside the lesions. For each lesion, the largest square area that could be completely 
enclosed in the lesion was selected, and an equivalent area of healthy skin outside the lesion 
was selected for comparison within the same image frame. The mean values of the p- 
parameters in Fig. 4 below were calculated in a similar way, by selecting a square area within 
each pathology that coincided with the region of higher p- parameter values. 
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Fig. 2. The mean optical density spectra of 17 examined melanomas (empty circles), 65 
common nevi (filled circles) and 82 healthy skin cases (dark filled circles). 

3. Results and discussion 

To illustrate the proposed differentiation parameter, Fig. 3 shows ordinary color RGB and 
parametric images of a melanoma and a nevus (a) and a pigmented nevus alone (b). The 

a) 
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Fig. 3. RGB images (left) and parametric p-images (right) of melanoma (a) and pigmented 
nevus (b). The p-value of healthy skin is set to 0. 
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average /?-value of healthy skin in these images is set to 0, in order to make the results 
comparable regardless of the individual skin type or color. The nevi regions in both cases are 
darker, while the melanoma region is lighter than the healthy skin background. 

Most of the parametric /?-maps of the examined lesions showed the same tendency — 
melanoma's areas (or parts of them) showed notably higher /7-values than the surrounding 
skin, and the nevus areas showed lower /^-values than the healthy skin, or, in some cases, 
similar p- values. Consequently, the criterion p > p 0 , where p 0 is related to the surrounding 
healthy skin, may be regarded as indication to melanoma. 

Figure 4 shows the scatter plot of the mean /^-values (averaged over the selected pathology 
areas) for all examined cases of melanoma and common nevus, when p 0 is set to 0. 

5 i 



a 
■ 



• melanoma 

o common nevus 



0 10 20 30 40 50 60 
case number 

Fig. 4. Scatter plot of the normed p- parameter for all melanomas (dark circles) and nevi (light 
circles). 

As the result of our diagnostic test, sensitivity, specificity, positive predictive value and 
negative predictive value were calculated. After analysis of the 17 available melanoma multi- 
spectral image sets, 16 showed clear response accordingly to the above-mentioned criterion p 
(true positives or TP), and only one case - melanoma with ulceration - did not support the 
criterion (false negative or FN). After analysis of 65 nevi, 58 of them showed lower p- 
parameter values than the surrounding skin, (true negatives or TN), and 7 cases did not meet 
the criterion (false positive or FP). Sensitivity of this test is TP/(TP + FN) = 94%, specificity 
is TN/(FP + TN) = 89%, positive predictive value (PPV) = TP/(TP + FP) = 70% and negative 
predictive value (NPV) = TN/(FN + TN) = 98%. 

4. Conclusions 

The reported results show that sensitivity and specificity of the proposed approach are high 
(94% and 89%, respectively). However, for clinical implementation it is necessary to verify 
this algorithm on a group of dysplastic nevi (the risk group of melanoma). But, presumably, 
there is a potential for clinical implementation, e.g. as additional test in cases when ABCD 
rule describing the four criteria of melanoma is not easy to apply, or in situations that cast 
doubts. Further studies would promote better understanding of the photophysiological 
processes in human skin and its lesions. 

The image processing technique can be further improved. For instance, a simpler form of 
Eq. (2) can be used: 



7(540) 



7(650) -7(950) 



(4) 
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wherC 7 o (65Q) - 7 o (95Q) . Besides, healthy skl n also can be taken as a reference, and then 
/ 0 (540) 

I 0 (540), Io(650) and Io(950) will be reflection intensities from normal skin. 

Improved software could replace the visual assessment (Figs. 2, 3) by automated selection 
of skin lesions in the future. 
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